First Principles Electronic Structure of Mn doped GaAs, GaP, 

and GaN semiconductors 

T. C. Schulthess 1 , W. M. Temmerman 2 , Z. Szotek 2 , A. Svane 3 , L. Petit 1 
Computer Science and Mathematics Division and Center for Nanophase Materials Sciences, 
Oak Ridge National Laboratory, Oak Ridge, Tennesee 37831-6164, USA 
2 Daresbury Laboratory, Daresbury, Warrington WA4 4 AD, UK 
^Department of Physics and Astronomy, 
University of Aarhus, DK-8000 Aarhus C, Denmark 
(Dated: May 13, 2010) 



1 



Abstract 

We present first-principles electronic structure calculations of Mn doped III-V semiconductors 
based on the local spin-density approximation (LSDA) as well as the self-interaction corrected local 
spin density method (SIC-LSD). We find that it is crucial to use a self-interaction free approach to 
properly describe the electronic ground state. The SIC-LSD calculations predict the proper elec- 
tronic ground state configuration for Mn in GaAs, GaP, and GaN. Excellent quantitative agreement 
with experiment is found for magnetic moment and p-d exchange in (GaMn)As. These results allow 
us to validate commonly used models for magnetic semiconductors. Furthermore, we discuss the 
delicate problem of extracting binding energies of localized levels from density functional theory 
calculations. We propose three approaches to take into account final state effects to estimate the 
binding energies of the Mn-d levels in GaAs. We find good agreement between computed values 
and estimates from photoemisison experiments. 

PACS numbers: 75.30.Et, 75.30.Hx, 75.50.Pp 
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I. INTRODUCTION 



Semiconductor spintronics aims to exploit both the spin and charge of electrons in new 
generations of fast, low dissipation, non- volatile integrated information storage and pro- 
cessing devices. The Mn doped Ga-V semiconductors are amongst the most interesting 
materials for applications in such new devices. In particular, (GaMn)As has been estab- 
lished as a well-behaved mean field ferromagnet with the Curie temperature T c linearly 
dependent on the concentration of the substitutional Mn,— a magnetic moment per Mn close 
to its free ion value,— and a possible low concentration of carriers further promoting T c .— 
Recent advances in thin-film growth of III-V semiconductors doped with Mn have led to 
synthesis of (GaMn)As dilute magnetic semiconductor (DMS) with Curie temperatures of 
the order of 173 K— . Also, recent experiments by Edmonds et al.— indicate a carrier-induced 
nature of the ferromagnetic exchange, but a small, finite, density of unoccupied Mn d states 
is found close to the Fermi level, reflecting hybridization with the host valence bands. Burch 
et al.,— on the other hand, claim observing impurity band conduction in Gai-^Mn^As, with 
large effective masses of the carriers, so far not confirmed by other experiments. 

The most interesting characteristic of DMS is the carrier induced nature of the magnetic 
coupling, which has two important practical implications. First, the carriers are polarized 
and DMS can serve as efficient sources for spin injection, owing to the fact that they are 
structurally compatible with semiconductors used in devices, which alleviates the problems 
of interfacial disorder that prohibits efficient spin-injection from traditional ferromagnets 
into semiconductors. Second, because the Curie temperature is correlated with the car- 
rier concentration, the magnetic order can be manipulated with voltage^. Therefore, efforts 
aimed at increasing the Curie temperature of magnetic semiconductors have to be concerned 
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with the nature of the magnetic exchange coupling in order not to loose these main advan- 
tages of carrier induced magnetism. However, in order that these materials be relevant for 
practical spin- and magneto-electronics^ applications, their Curie temperatures have to be 
raised above room temperature. 

The successful description of magnetic properties in Gai-^Mn^As motivated Dietl et 
al.— to use the Zener mode l 10 ! 11 description to predict Curie temperatures of various Mn 
doped group IV, III-V, and II- VI semiconductors. In particular, their prediction of a high 
Curie temperature of Mn doped GaN inspired many groups to synthesize this system, and 
several reports of Curie temperatures well in excess of room temperature now exist in the 
literature^. However, despite important advances made for the Gai-^Mn^As system, the 
microscopic nature of the electronic structure and magnetic exchange of Ga-V systems, and 
especially (GaMn)N, is far from understood. In particular, it still needs to be established 
whether the magnetism in these materials complies with the description provided by the 
Zener model. More fundamentally, it is not clear to what extent the Kondo-like Hamiltonian, 
describing Mn spins of 5=5/2 interacting with free carriers, is a justified starting point 
for describing systems other than Mn doped GaAs. For example, electron spin resonance 
measurements, which indeed support the picture of divalent Mn in GaAs^, clearly favour 
the trivalent d 4 configuration for Mn impurities in GaP— . For GaN it has also been reported 
that Mn is in a divalent state when electrons are dopecU^, but in a trivalent state when holes 
are doped to the system— Hwang et al.,— using photoemission and soft X-ray absorption 
spectroscopy, confirm that in the n-type doped GaN the Mn state is divalent, while for the 
non-doped one it is trivalent. 

Most model descriptions of Mn doped GaAs assume the Mn impurity to have a local- 
ized moment of S — 5/2, formed by five occupied atomic- like ci-orbitals (d 5 ) that interact 
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weakly with the host valence band through level repulsion, leading to the simple (d 5 + h) 
picture where the Mn spin couples antiferromagnetically to a polarized hole, h. The result- 
ing Kondo-like Hamiltonian appears well justified with parameters that can be established 
experimentally. However, its solution is still under debate. Dietl et al.^£, as well as Mac 
Donald and co-workers 1 - use the Zene r 10 ' 11 description, in which it is assumed that the 
ferromagnetic interaction between Mn ions is mediated by the induced spin-density that is 
anti-aligned with the Mn moments, and the resulting prediction of the Curie temperature, 
T c , appears in agreement with experiment for Gai-^Mn^As, with x up to 0.09 4 . 

The key question that needs to be addressed from an electronic structure point of view 
is whether at concentrations of several atomic percent of substitutional Mn, the acceptor 
level forms an impurity band broad enough to merge with the valence band. The electronic 
structure of Mn in Ga-V semiconductor hosts is best studied with first principles elec- 
tronic structure calculations. Since the Local (Spin-) Density Approximation 2 ^ (L{S}DA) 
to Density Functional Theor y 21 ' 22 (DFT) has proven very successful in predicting ground 
state properties in semiconductors 2 ^ as well as transition metals and their alloys 2 ^, several 
group o 25 ! 26 ' 27 ! 28 ' 29 ! 30 have performed LSDA calculations of the electronic structure and mag- 
netism of various transition metal doped semiconductor systems soon after the discovery of 
ferromagnetism in Mn doped GaAs. These studies all agree in one important point 2 ^: the 
Mn ci-orbitals are not atomic-like but hybridize rather strongly with the host valence band 
(As-p in the case of Gai_ x Mn x As). In this, the LSDA picture differs substantially from the 
one of model calculations. We have recently shown^i, that first principles calculations based 
on the self-interaction corrected local spin density (SIC-LSD) method, lead to a picture that 
is more consistent with what is expected for a system with strongly correlated electrons. The 
purpose of the present paper is to make a more in-depth comparison between the electronic 
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structure results using different functional and comparing the first principles predictions 
with experiments. We begin with a discussion of the SIC-LSD method, how it is applied in 
Mn doped III-V systems and how it compares to LDA and other methods used to "fix" LDA 
with respect to electron correlations. Then we compare the results of the SIC-LSD calcula- 
tions for Mn doped GaAs with experiment and previously known LDA results. Finally, by 
studying the changes in electronic configuration and trends for different III-V systems, we 
shed some light on the intricacy of the electronic structure of Mn doped GaAs. 

II. SPURIOUS SELF-INTERACTIONS AND THEIR REMOVAL 

The local (spin) density approximation to density functional theory forms the founda- 
tion of the conventional band theory. It describes electron correlations at the level of the 
homogeneous interacting electron gas and has been very successful in predicting electronic 
properties of many materials in terms of their ground state charge densit y 32 > 33 For semi- 
conductors, this is particularly true for the prediction of structural properties and alloying 
behavior. A common criticism of LDA and its generalization, the generalized gradient ap- 
proximation (GGA), namely the systematic tendency to underestimate the band gap, is 
strictly speaking not a shortcoming, since DFT in its present use is a theory of the ground 
state and hence not meant to predict quasi particle excitation spectra. Nevertheless, LDA 
and GGA based band structures are often used to describe the valence and conduction bands 
in semiconductors and numerous methods have been applied to correct the magnitude of 
the band gap.— 

In solids with localized d- and/or /-electrons, such as the transition metal monoxides, the 
cuprate high temperature superconductors, rare earths, actinides, and, as we will see below, 
dilute magnetic semiconductors, LDA notoriously fails to describe correctly the electronic 
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and magnetic ground state properties. This can be understood as follows: in Kohn-Sham 
based DFT, the energy functional of the electron density is written as, 

E[n\ = T s [n] + J V ion (r)n(r)dr + J[n] + E xc [n], (1) 

where T s is the single particle kinetic energy, V lon the ionic potential, J the classical Coulomb 
interaction energy - or Hartree energy - of the charge density, and E xc is the exchange 
and correlation energy that contains the exchange term and all the non-classical electron 
correlations. When the electron density is decomposed into orbital densities, n = Yli n ij 
it is straightforward to demonstrate that the Hartree term contains a contribution, J[ni\, 
of an orbital interacting with itself. This self-interaction term is cancelled exactly by the 
self-exchange contribution to E xc . In the LDA and GGA, E xc is approximated and the 
self-interactions are not cancelled anymore. One speaks of spurious self-interactions that are 
introduced by the approximations in the LDA and GGA functionals. 35 

The spurious self-interactions are negligible for extended orbitals such as the s and p 
bands in semiconductors or d bands in transition metals. They are, however, substantial 
whenever electrons occupy localized orbitals such as the 3d orbitals of transition metal atoms 
in oxides or of transition metal impurities in semiconductors. In these cases, the spurious 
self-interactions push the localized orbitals into the valence band usually resulting in too 
strong a hybridization with the other valence electrons. 

This problem was recognized already more than two decades ago and a remedy was 
proposed by Perdew and Zunger 36 to simply subtract the spurious self-interactions from the 
LSDA functional. Their self-interaction corrected (SIC) local spin density (LSD) functional 
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takes the form 



occ. 

£SIC = £LSDA_^SIC (2) 

i 

where the sum runs over occupied and localized orbitals with non- vanishing self-interaction 
corrections 

5f c = JN + ^ c SDA N- (3) 

When applied to atoms, the most extreme case where all electrons occupy localized 
orbitals, the SIC-LSD functional drastically improves the description of the electronic 
structure*^. In solids, where not all electrons occupy localized orbitals, one is faced with 
the task of minimizing the orbital dependent SIC-LSD functional (J2J). Additionally, since 
the LSDA exchange correlation functional depends non-linearly on the density, the self- 
interaction corrections (JHJ) and hence the SIC-LSD functional (J2J) are not invariant under 
unitary transformations of the basis and one is thus faced with a daunting functional mini- 
mization problem. 

Since the main effect of the self- interact ion correction is to reduce the hybridization of lo- 
calized electrons with the valence band, the technical difficulties of minimizing the SIC-LSD 
functional in solids can often be circumvented by introducing an empirical Coulomb inter- 
action parameter U on the orbitals that are meant to be localized. The original derivation 
of the LDA+U approach^ seems to have been based on the conjecture that LDA can be 
viewed as a homogeneous solution of the Hartree-Fock equations with equal, averaged, occu- 
pations of localized d- and/or /- orbitals in a solid. Therefore, as such, it can be modified to 
take into account the on-site Coulomb interaction, U, for those orbitals to provide a better 
description of their localization. The on-site Hubbard U is usually treated as an adjustable 
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parameter, and chosen to optimize agreement with experiment. In many transition metal 
oxides, the approach is successful because the results are not very sensitive to the precise 
value of U . However, in Mn doped III-V systems, the opposite seems to be the case, 39 
magnetic exchange and moment depend sensitively on the value of U and the method looses 
its predictive power. 

For this reason we are using an approach 4 ? with which the full SIC-LSD functional is 
minimized with respect to the orbital decomposed charge density, giving rise to a generalized 
eigenvalue problem with an orbital dependent potential - as the self-interaction corrections 
are only non-zero for localized electrons^, the localized and delocalized electrons experience 
different potentials. The latter move in the LSD potentials, defined by the ground state 
charge density of all occupied states, while the former experience a potential from which 
the self-interaction term has been subtracted. Hence, in this formulation one distinguishes 
between localized and itinerant states and it is possible to study different nominal valences 
for those elements in the solid that contain localized d— and/or /-electrons. To determine 
the ground state energy and valence, one minimizes the SIC-LSD functional with respect to 
these electronic configurations (different distributions of localized and itinerant states). The 
resulting SIC-LSD method is a first principles theory for the ground state with no adjustable 
parameters. Finally, it is important to note that the SIC-LSD functional subsumes the 
LSDA, that is, when all electrons (besides the core electrons) are itinerant, the configuration 
in which no orbitals self-interaction correct will have the lowest energy and the solutions of 
the SIC-LSD functional will be identical to that of the LSDA. 
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III. ELECTRONIC STRUCTURE OF GAMNAS 



In this study we model the idealized 4 ^ Mn doped semiconductor with a supercell approach 
with a unit cell of up to 64 atoms, in which one Ga atom has been replaced with a Mn 
impurity. We limit ourselves to the zinc-blende structure with experimental lattice constants 
of the host semiconductor and neglect structural relaxations. In all calculations, the 4s and 
4p electrons of the semiconductor hosts are itinerant, forming the valence band - applying 
self-interaction corrections to these electrons raises the energy. Self-interaction corrections 
are applied to states that originate from the Mn-3<i electrons, that is, we minimize (J2J) 
over orbital decomposed densities and the different combinations of applying self-interaction 
corrections to the Mn-3<i levels. Of all possible combinations, the following two scenarios 
always have the lowest energy: (1) when self-interaction corrections are applied to all five 
majority Mn-3<i levels - this is the 5 = 5/2 configuration which we denote by Mn(<i 5 ) or 
Mn 2+ ; (2) when self-interaction corrections are applied to all but one of the majority Mn 
t2 g levels - this configuration has S = 2 and we will call it Mn(d A ) or Mn 3+ . 

A. Ground state, magnetic moment, and exchange 

The relative energies of the three relevant configurations for one Mn impurity in a 64 
atom unit cell of GaAs are E(d A ) - E(d 5 ) = 0.18 eV and E(LSDA) - E(d 5 ) = 1.99 eV. The 
lowest energy state (ground state) of this system is reached when all five majority d electrons 
localize on the Mn atoms - the same is true for higher Mn concentrations (smaller unit 
cells). The ground state predicted by these calculations agrees with electronic configuration 
inferred from electron spin resonance measurements^ . Both the d 5 and the d A scenarios are 
favourable when compared to the LSDA solution (for which none of the five Mn d orbitals 
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localize). The magnetic moment on the atomic sphere surrounding the Mn atoms for the 
d 5 , d 4 , and LSDA scenarios is 4.50/ie, 4.07/ib, and 3.80/ib respectively. The value for the 
ground state is in excellent agreement with the local Mn moment measured experimentally 
using X-ray magnetic circular dichroism (XMCD). 2 

A detailed characterization of the electronic state of Mn impurities in GaAs and the 
magnetic coupling can best be seen from the density of states (DOS), which is plotted in 
figure n In the d 5 ground state scenario, all five majority <i-states form a localized impurity 
band below the host valence band with virtually no hybridization. However, because of 
repulsion between the majority Mn-d levels and the As-p bands, the top of the host valence 
band is spin split, leading to a spin polarized hole. In this scenario, the hole has primarily 
As-p character. This is a first principles electronic structure description of what is called 
the (d 5 + h) state in the models, where the Mn impurity with S = 5/2 is surrounded by a 
hydrogenic hole of opposite spin polarization. The spin polarization of the hole is, however, 
not 100%, since, as can be seen from the plot, the Fermi level crosses both the majority 
and the minority bands. In other words, the spin splitting A of the valence band is smaller 
than the difference, -Evbm — -^F; between the valence band maximum and the Fermi level. 
That the hole is not strongly coupled to the Mn impurity, can also be inferred from its 
spacial extent: in our calculations the induced hole in the valence band extends beyond 
the size of the unit cell. Hence, contrary to what has been assumed in some models of 
Mn doped GaAs, the coupling between the Mn S = 5/2 spin and the host valence band is 
not strong, and the Zener model based approach 9,1 ^ to estimating the Curie temperature is 
justified. From the splitting of the valence band maximum we can estimate the p—d exchange 
coupling (see ref.— for details), the predicted values are given in Tabled The agreement with 
the experimental value extracted from spin polarized photoemission dat a 42 ' 43 is very good. 
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The (d 5 + h) picture for Mn doped GaAs seems to be well justified by the first principles 
calculations. However, this applies only in the concentration range below 10 at. % Mn. 
In calculations with 8 -atom unit cells, corresponding to Mn concentrations of 25 at. %, 
respectively, the valence band is strongly perturbed by the Mn alloying effect and the p — d 
exchange parameter extracted from the valence band splitting is concentration dependent. 
The simple Zener model is probably not valid at higher Mn concentrations and one should 
not simply extrapolate the Curie temperature predictions of this model to concentrations 
higher than 10 at. % Mn. 

The effects of the spurious self-interactions in the LSDA are most clear in the density of 
states [TJ Since localizing the Mn-d orbitals is energetically unfavourable due to the unphys- 
ical self-interactions, the d levels are pushed up into the valence band and hence become 
strongly hybridized. As a result, the local Mn moment (3.8 hb for LSDA) is underestimated. 
But more importantly, the nature of the states at the Fermi level is very different leading to 
a different mechanism for the magnetic coupling. The p — d exchange coupling is strongly 
overstated, the splitting of the top of the valence band being a factor of two larger than 
experiment (Table 

B. Localized levels and photoemission 

Although DFT is strictly speaking a theory of the ground state from which spectroscopic 
information is not easily extracted, the LDA based band-structure is often compared to pho- 
toemission experiments. This is because the effective Kohn-Sham potentials can be viewed 
as an energy independent self-energy and hence the Kohn-Sham energy bands correspond 
to the mean field approximation for the spectral function. In the SIC-LSD, this argument 
only applies to the itinerant orbitals that are not self-interaction corrected. The localized 
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FIG. 1: The spin-resolved total DOS (solid red line) and Mn partial DOS (dotted blue line) of 
(Gai-^Mn^As with x = 6.25% for the SIC-LSD ground state, d 5 , scenario (upper panel) as well 
as the LSDA (lower panel). 

states that have been self-interaction corrected see a different potential, and the solution 
(or the solutions) to the generalized SIC-LSD eigenvalue problem, which is different from 
the solution to the Kohn-Sham equations in the LDA, no-longer correspond to a mean field 
approximation of the spectral function. To extract spectroscopic information for the local- 

13 



TABLE I: Spin splitting at the top of the valence band, A and the p — d exchange coupling 
parameter J p d and the more commonly quoted Nq/3 parameter (Nq is the number of unit cells in 
the normalization volume and = —J p d). We are comparing the ground states of the SIC-LSD 
and the LSDA functionals with experiment. 





X 


A 


Jpd 


N p 






(eV) 


(meV nm 3 ) 


(eV) 


Mn((i 5 ) 


0.25 


0.69 


27 


-0.60 




1/16 


0.43 


68 


-1.51 




1/32 


0.21 


67 


-1.49 


LSDA 


1/16 


0.71 


114 


-2.5 




1/32 


0.45 


145 


-3.2 


Experiment 






54±9 


-1.2±0.2 



ized Mn 3d levels from the SIC-LSD calculations, we have to take a different route. We 
consider three approaches: 

1. The canonical approach to compute photoemission energies of localized states in 
the literature is a Ascf calculation^, in which the total energy differences between two 
configurations of Mn assumed to represent the initial and final states of a photoemission 
experiment are compared. This is expressed by the formula 

A^scf = E(d\ N pd + 1) - E(d 5 , N pd ), (4) 

meaning that one removes an electron from a Mn d level and introduces it at the top of 
the valence band, thus increasing the total number of valence electrons in the supercell, 
N p d, by 1. Here E represents the total energies of the respective d A and d 5 configurations. 
For E(d A , N p d + 1) the electronic configuration is constrained such that the fifth d level is 
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TABLE II: Positions of the Mn 3d level (in eV) as they would appear in photoemission experiments 
computed with the methods discussed in the text. 



A (l) 

^SCF 


A (2) 
^SCF 


OEP 




experiment 


3.23 


3.16 


3.7 


5.1 


4.2 



unoccupied. Of course, a problem one faces with this kind of estimate is that the correct 
final state might not be the ground state of the system with one electron removed. For 
example, if a final state of d is assumed for Mn in GaAs, and the constraint on the fifth 
d level is dropped during iterations to self-consistency it will move below the valence band 
maximum (VBM), and the level starts to get filled (until eventually the filling balances the 
position of the level at the Fermi energy). Hence, in the end the final state corresponds to 
some intermediate Mn configuration between d 5 and constrained d A . The effect is partially 
physical, since the screening processes around a hole created by photoexcitation include Mn 
d screening, however not by the state being kicked out but by all the minority d electrons 
which become attracted to the <i-hole. An alternative interpretation of the photoemission 
experiment is given by the formula 

A^scf = E(d\ N p d) - E(d\ N pd - 1). (5) 

where one simply compares two final state energies: the energy of the system with one Mn 
d electron removed compared to the energy of the system with one electron removed at the 
Fermi level. This alternative way of interpreting the photoemission spectrum leads to the 
numerical result that is very similar to Ag^ F , as can be seen from Table ILT1 

2. In the earlier section we have described the SIC-LSD method as an orbital depen- 
dent density functional theory. Formally, however, SIC-LSD may be viewed as a standard 
density functional theory, implying that the SIC-LSD energy functional can be represented 
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as a functional of the total charge density alone and minimized with respect to it. This 
means that there exists an effective Kohn-Sham equation with an effective potential, which 
is common to all Kohn-Sham states, is self-interaction free, and depends only on the to- 
tal charge density. Here the situation is completely analogous to the optimized effective 
potential (OEP) introduced in connection with the Hatre-Fock approximation , 45 ! 46 ! 47 ' 48 for 
which case the Kohn-Sham eigenvalues are often compared to quasiparticle energies, with a 
considerable improvement over the Hartree-Fock eigenenergies. 

Adopting the OEP philosophy, one can search for the effective potential, which reproduces 
the SIC spin-density, and with such a potential derive the density of states of Mn in GaAs. 
Since the only significant difference between the LSD and SIC-LSD charge densities for Mn 
in semiconducting hosts is the Mn spin moment, we can constrain the search by looking for 
that particular potential shift on the Mn site, which will reproduce the self-consistent spin 
moment of Mn in the SIC-LSD calculations. Hence, as in the standard OEP approach, the 
derived eigenvalues reflect the self-interaction correction in that the Mn majority d-states lie 
lower in energy than in the LSD case, however not as low as when calculated directly from 
the SIC eigenvalues. The OEP-like SIC-derived density of states is shown in figure 121 One 
can see nearly perfect agreement with the experimentally determined position of the Mn-d 
states ( Table ITT]). 

3. As a 'quick and dirty' way to obtain electron removal energies just from the self- 
consistent ground state calculation one may apply a transition state approximation, accord- 
ing to which the Mn removal energy is defined as the average of the SIC-LSD and LSD d 
state positions: 

e TS = \{< /I^lsd + Vsic\f > + < /I^lsdI/ >)• (6) 
In effect, the SIC potential is only counted with half of its strength in the transition state 
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FIG. 2: Density of states from SIC-based OEP calculations, showing both the total (solid red line) 
and Mn d (dotted blue line) contributions. 

approximation to the removal energy. By evaluating Hlsd in the initial state, i.e. without 
the hole in the d shell, we avoid the aforementioned effect of the <i-hole pulling the d levels 
down. This approximation has been discussed for rare-earth imp urities in GaAs and GaN.— 



The transition state philosophy was also implemented in Ref. 



5fj, albeit in a different manner, 



by invoking the averaging factor of \ already in the total energy functional, while we do it 
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here only for the removal energy, Eq. (jHJ), after self-consistency. 

In Table HT1 we compare the three methods just described to estimate the positions of the 
d levels in photoemission experiment of Mn doped GaAs. The agreement with experiment 
is satisfactory, but certainly not as perfect as the quantitative agreement we find for the 
magnetic moment (discussed above). This should, however, come as no surprise, since DFT 
is a theory for the ground state and the magnetic moments that are determined from the 
spin densities are rigorously founded in spin depended DFT. Similar results for the d levels 
have been achieved with the LDA+U method by Shick et al.,— assuming U of 4 eV. Note, 
however, that in that study they report a large sensitivity of the Mn d binding energies to 
the magnitude of U, which is not the case in the present work, as there is no adjustable 
parameter in the SIC-LSD method. 

IV. MATERIALS TRENDS: MN IN GAAS, GAP, AND GAN 

The energy differences, E(d 4 ) - E(d 5 ), between the Mn(ci 4 ) and Mn(c/ 5 ) scenarios in a 64 
atom unit cell of GaAs, GaP, and GaN are, respectively, 0.18 eV, -0.05 eV, and -1.28 eV. 
Hence, the ground state configuration of Mn changes from d 5 in GaAs to d A in GaP and 
GaN. While the energy difference between the configurations is very small for Mn doped 
GaP, the predicted d A ground state seems to be in agreement with electron spin resonance 
experiments^. This change in valence of the Mn impurity can best be understood in terms of 
a localization/delocalization transition of the fifth Mn d orbital in the majority channel. In 
figure|3]we compare the density of states of the two scenarios of Mn in all three Ga-V systems 
studied here. In GaN, a charge transfer insulator, the fifth Mn d level is unoccupied, forming 
a deep impurity level in the middle of the band gap, and there is virtually no hybridization 
with the valence band. In GaP, the lattice constant is increased, the band gap reduced, the 
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valence band moves up and closer to the Mn d impurity level, and hybridization between the 
two is already substantial. As the hybridization in the d 4 scenario further increases in GaAs, 
the occupation of the fifth Mn d level has increased to the point, where, due to the larger 
lattice constant, it is now energetically more favorable to localize, leading to a d 5 ground 
state configuration for the Mn impurity and a hole with As p character. 

Since the SIC-LSD method is a static theory for the ground state, our present calculations 
do not account for fluctuations that are almost certainly playing an important role whenever 
two or more solutions are close in energy, i.e. close to the transition between localized and 
delocalized states. Quantum fluctuations as well as fluctuations at finite temperature may 
well lead to some admixture of Mn d A and Mn d 5 in GaAs and GaP. Furthermore, co-doping 
with additional donors or acceptors will likely also have an effect on the localization nature 
of the fifth Mn d orbital.— Additionally, since interstitial Mn impurities, which we did not 
consider here, give rise to a localized impurity band near the Fermi level^i, the measured 
electronic structure in Mn doped GaAs can be affected by many factors that are sample 
dependent. 

V. SUMMARY 

As we have seen, the use of standard band structure techniques, such as the LSDA, does 
not lead to satisfactory description of the electronic structure of Mn doped semiconductors. 
This is because spurious self-interactions introduced by the approximation to the exchange 
correlation functional lead to an incorrect description of the Mn d orbitals as well as of the 
hole mediated magnetic exchange. Hence, to properly describe magnetic semiconductors, 
methods have to be used that do not suffer from spurious self-interactions. 

Here, we have applied the SIC-LSD method to study the electronic structure and magnetic 
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FIG. 3: The spin-resolved total DOS (solid red line and bottom scale) and Mn-DOS (dotted blue 
line and top scale) from left to right, for Mn-d 4 impurity in GaN, GaP, and GaAs, as well as Mn-ci 5 
impurity in GaAs. Results are shown for 64 atom supercell (3.125% Mn). 

properties of Mn doped III-V semiconductors. For (GaMn)As our calculations predict the 
correct electronic ground state configuration for Mn, with its magnetic moment of 4.5 
and p-d exchange in excellent agreement with experiment. In agreement with previously 
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assumed models, we find that the Mn spin is weakly antiferromagnetically coupled to the 
hole that mediates the ferromagnetic exchange. By taking into account screening/relaxation 
effects, we have obtained a very good agreement of the calculated binding energies for the 
Mn d states with photoemission experiments. Our calculations also describe correctly the 
trends of the electronic ground state for Mn in GaN, GaP, and GaAs. In particular we 
predict the localization/delocalization transition from Mn (d 5 ) to Mn (d 4 ), when moving 
from GaAs to GaP and GaN. As the Mn(rf 5 ) and Mn(d 4 ) configuration are energetically 
very close in GaP and GaAs, it is likely that effects of co-doping, additional impurities, and 
thermal fluctuations lead to mixed valence states for the measured electronic ground state 
configuration of Mn that might be sample dependent. 
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